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For Anderson localization on the Cayley tree, we study the statistics of various observables as 
a function of the disorder strength W and the number A'^ of generations. We first consider the 
Landauer transmission Tjv. In the localized phase, its logarithm follows the tra veling wave form 
QQ ' In Tjv — In Tjv + Int* where (i) the disorder- averaged value moves linearly In(rjv) — a^nd the 

■ localization length diverges as ^loc ~ {W — Wc)'"'"" with uioc ~ 1 (ii) the variable t* is a fixed random 

O ■ variable with a power-law tail P*{t*) ~ for large t* with < l3{W) < 1/2, so that 

all integer moments of Tjv are governed by rare events. In the delocalized phase, the transmission 
^ . Tjv remain s a finit e random variable as A — > oo, and we measure near criticality the essential 

Q ' singularity In(roo) ~ —\Wc — Wl'"^^ with kt ~ 0.25. We then consider the statistical properties 

of normalized eigenstates — !> in particular the entropy S = — |V'(2^)P In 

and the Inverse Participation Ratios (I.P.R.) Iq = \iIj(x)\'^'^ . In the localized phase, the typical 
' entropy diverges as Styp ~ {W — Wc)~"^ with vs ~ 1.5, whereas it grows linearly Styp{N) ~ A^ in 

the delocalized phase. Finally for the I.P.R. , we explain how closely related variables propagate as 
traveling waves in the delocalized phase. In conclusion, both the localized phase and the delocalized 
phase are characterized by the traveling wave propagation of some probability distributions, and 
the Anderson localization/delocalization transition then corresponds to a traveling/non-traveling 
critical point. Moreover, our results point towards the existence of several length scales that diverge 
C/3 . with different exponents v at criticality. 



. I. INTRODUCTION 

S: 

' ■ Since its discovery fifty years ago [1] Anderson localization has remained a very active field of research (see for 
y instance the reviews [2-7]). According to the scaling theory [8], there is no delocalized phase in dimensions d = 1,2, 
Q vifhereas there exists a localization/delocalization at finite disorder in dimension d > 2. To get some insight into this 
^ ] type of transition, it is natural to consider Anderson localization on the Cayley tree (or Bethe lattice) which is expected 
to represent some mean-field limit. The tight-binding Anderson model on the Cayley tree has been thus studied by 
^ [ various authors. In [9], the recursion equation for the self-energy has been studied to establish the mobility edge as 
^ the limit of stability of localized states. In [10], the exponent v governing the divergence of the localization length 
' was shown to be = 1. In [11], the supersymmetric formalism has been used to predict the critical behavior of some 
P\J , disorder-averaged observables (see also [12-15] where similar results have been obtained for the case where Anderson 
■ tight-binding model is replaced by a non-linear cr-model). In [16, 17], the problem was reformulated in terms of 
'""^ , recursions on Riccati variables to obtain weak-disorder expansions. Other studies have focused on random-scattering 
I ■ models on the Cayley tree [18-20]. More recently, the interest in Anderson localization on the Cayley tree has been 
T-H ' revived by the question of many-body localization [21], because the geometry of the Fock space of many-body states 
was argued to be similar to a Cayley tree [22-25]. 

In the present paper, we consider the Anderson tight binding model on the Cayley tree already studied in [9- 
^ , 11, 16, 17] and we study numerically the statistical properties of Landauer transmission and eigenstates as a function 
of the disorder strength W and of the number N of generations. We find that several probability distributions 
^ . propagate as traveling waves with a fixed shape, as the number N of tree generations grows. The fact that traveling 
• waves appear in disordered models defined on trees has been discovered by Derrida and Spohn [26] on the specific 
example of the directed polymer in a random medium and was then found in various models [27]. For the case of 
Anderson localization on trees, the analysis of [9] concerning the distribution of the self-energy in the localized phase is 
actually a 'traveling wave analysis' (see Appendix A for more details), although it is not explicitly mentioned in these 
terms in [9]. The fact that traveling waves occur has been explicitly seen in the numerical study of the transmission 
distribution for a random-scattering model on the Cayley tree (see Fig. 8b of [20]) and has been found within the 
supersymmetric formalism [12, 14, 15]. In the field of traveling waves and front propagation (see the reviews [28, 29]), 
there exists an essential separation between two classes : in 'pulled fronts', the velocity is determined by the form 
of the tail of the front and thus by the appropriate linearized equation in the tail region, whereas in 'pushed fronts', 
the velocity is determined by the bulk properties and thus by the non-linear dynamics in the bulk region. It turns 
out that for disordered systems defined on trees, the traveling waves that appear usually correspond to cases where 
it is the tail of the probability distribution that determines the velocity [26]. In particular for Anderson localization 
on the Cayley tree, the traveling wave propagation is also determined by the tails [9]. From this traveling wave point 
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of view, the localization/delocalization Anderson transition thus represents a traveling/non-traveling critical point. 
Such a traveling/non-traveling phase transition for a branching random walk in the presence of a moving absorbing 
wall has been studied recently in [30] and we find here very similar behaviors in the critical region. 

The paper is organized as follows. In Section II, we describe the statistical properties of Landauer transmission. In 
Section III, we discuss the statistical properties of eigenstates, as measured by the entropy and Inverse Participation 
Ratios (I.P.R.). We summarize our conclusions in Section IV. In Appendix A, we translate the analysis of Ref. [9] 
concerning the distribution of the self-energy in the localized phase into a traveling-wave tail analysis for the Landauer 
transmission discussed in section II. In Appendix B, we explain how such a similar tail analysis can be performed 
in the dclocalizcd phase for auxiliary variables that arc closely related to Inverse Participation Ratios. Finally, in 
Appendix C, we recall the results of Ref. [30] concerning the traveling/non-traveling phase transition for a branching 
random walk in the presence of a moving absorbing wall, since these results are used as a comparison in the text to 
understand the finite-size scaling properties in the critical region. 



II. STATISTICAL PROPERTIES OF THE LANDAUER TRANSMISSION 



In quantum coherent problems, the most appropriate characterisation of transport properties consists in defining a 
scattering problem where the disordered sample is linked to incoming wires and outgoing wires and in studying the 
reflexion and transmission coefficients. This scattering theory definition of transport, first introduced by Landauer 
[31], has been much used for one-dimensional systems [32-34] and has been generalized to higher dimensionalities 
and multi-probes measurements (see the review [35]). In dimension d = I, the transfer matrix formulation of the 
Schrodinger equation yields that the probability distribution of the Landauer transmission becomes asymptotically 
log-normal [32, 34] , i.e. one has 

lnTl'''' = -A + Li/2« (1) 

where ^joc represents the localization length and where u is a sample-dependent random variable of order 0(1) 
distributed with a Gaussian law. Although it is often assumed that this log-normal distribution persists in the 
localized phase in dimension d = 2,3, recent numerical studies [36] are in favor of the following scaling form for the 
logarithm of the conductance 

In5f = -A+L-W„ (2) 

with exponents of order Lo{d = 2) ~ 1/3 and u){d = 3) ~ 1/5 [36], whereas u){d = 1) = 1/2 from Eq. 1. For the Cayley 
tree that we consider in this paper, we will find below that the fluctuation exponent vanishes ivcayiey = 0, and we 
will discuss the probability distribution of the variable u. But let us first recall the appropriate scattering framework 
for the Cayley tree [17] . 



A. Reminder on the Miller-Derrida framework to compute the Landauer transmission [17] 

We consider the Anderson tight-binding model 

iI = ^ei|z><il-F ^ jiXil (3) 

i <i,j> 

where the hopping between nearest neighbors < i,j > is a constant V = 1 and where the on-site energies e, are 
independent random variables drawn from the flat distribution 

The parameter W thus represents the disorder strength. 

We consider the scattering geometry introduced in [17] and shown on Fig. 1 : the finite tree of branching ratio K 
is attached to one incoming wire at its root (generation n = 0) and to K"^^ outgoing wires at generation 2N. One is 
interested into the eigenstate \tp > that satisfies the Schrodinger equation 



H\iIj >= Ftp > 



(5) 
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iljj{n>i) = tje*''("-*^ 



FIG. 1; Scattering geometry of Ref. [17] : the disordered tree of branching ratio K = 2 starting at generation n — and 
ending at generation 2N (on the Figure 2N = 4) is attached to one incoming wire and to K^^ outgoing wires. In section II, 
we discuss the properties of the total transmission T = \tj\'^ = 1 — \r\^ where r is the reflexion amplitude of the incoming 
wire, and tj the transmission amplitudes of the outgoing wires. 



inside the disorder sample and in the wires where one requires the plane-wave forms 

^{n < 0) 



V'j(n > 2N) 



ik{n-2N) 



(6) 



These boundary conditions define the reflexion amplitude r of the incoming wire and the transmission amplitudes tj 
of the j = 1, 2, ..K^^ outgoing wires. To satisfy the Schrodinger Equation of Eq. 5 within the wires with the forms 
of Eq. 6, one has the following relation between the energy E and the wave vector k 

E = 2cosk (7) 

To simplify the discussion, we will focus in this paper on the case of zero-energy E — and wave- vector k = t:/2 

In this paper : E = and k = 7r/2 (8) 

because the zero-energy E = corresponds to the center of the band where the delocalizion first appears when the 
strength W of the disorder is decreased from the strong disorder localized phase. 

Inside the Cayley tree < n < 2N — 1, the Schrodinger Equation of Eq. 5 involves one ancestor denoted by 
anc{n,j) and K descendants denoted by desm{n,j) 



K 







-{■n,j)ip{n,j) +il){anc{n,j)) + ^ il){desm{n, j)) 



m=l 



whereas for the last generation 2N there is only one ancestor and one descendant (outgoing wire) 

= €{2N,j)ilj{2N,j)+il}{anc{2N,j))+ilj{2N + 
As explained in [17], it is convenient to introduce the Riccati variables 



(9) 



(10) 



(11) 



that represent the ratio of the wave function of two neighboring sites. On the outgoing wires, these Riccati variables 
are fixed by Eq. 6 to be 



R{2N + l,j) = 



ij{2N+l,j] 



= e-''' = -i 



(12) 
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The Schrodinger equation of Eq. 10 gives the first recursion 

Ri2N, j) = -ei2N, j) - j^^^j^^^^^^ = -e(2A^, j) - e'^ = -e(27V, j) -i (13) 
whereas the Schrodinger equation of Eq. 9 gives the recursion inside the tree for < n < 2N — 1 

«(.M) = -«(.M)-f J^^Sj;^ (14) 

On the other hand, the value of the Riccati variable for the origin of the tree n = is fixed by the incoming wire of 
Eq. 6 

^ ' V(0) l + r ^ ' 

and thus the reflexion coefficient r can be obtained via 

i?(0) - e'fc I - R{0) ^ ' 

from the i?(0) obtained via the recursion of Eq. 14. From the conservation of energy, the total transmission T is 
related to the reflexion coefficient |rp 

T^^li.f = 1-H2 (17) 

3 

As explained in [17], the criterion for the localization/delocalization phases is then the following : 

(a) if the Riccati variable i?(0) at the root of the tree converges towards a real random variable as A'' — > oo, the 
reflexion is total |r| = 1 and the transmission vanishes T = 0. 

(b) if the Riccati variable i?(0) keeps a finite negative imaginary part as A — > oo, the reflexion is only partial \r\ < 1 
and the transmission remains finite T=l — |rp>0 

We refer to [17] for the results of a weak disorder expansion within this framework, and for a numerical Monte-Carlo 
approach to determine the mobility edge in the plane {E, W). Here we will instead study the statistical properties of 
the transmission T at zero energy E = as, a. function of the disorder strength W (Eq 4) and of the number N of 
generations. But before discussing the disordered case, let us first describe the finite-size properties of the pure case. 



B. Example : transmission of finite pure trees 



In the pure case where all on-site energies vanish e{n,j) = 0, all branches are equivalent and there is no dependence 
on j. The recursions for Riccati variables of Eqs 12, 13 and 14 give 



Rp^rei^N+l) =-i 
Rpure{2N) = -i 

K 

Rpureij^ 

i.e. one obtains the simple alternation between odd and even generations 

K 

Rpure{2N) 

K 

RpureiP'N — 1) 

(19) 

-Rpt,re(l) = -Ki 

(20) 



Rpurein) =-- Ti-T-TT for < n < 27V - 1 (18) 



Rpurei2N - 1) = = -Ki 



Rpure{2N - 2) = J^^^f-^ = -i 
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The reflexion coeflicient r thus vanishes exactly for any even size (2A'') 

_ i + Rpure{0) 



and the transmission coeflicient in each branch j = 1, .., K"^^ reads 



= (21) 



2JV 



(Note that for pure trees of uneven size, the reflexion would not vanish. This is why in the disordered case, we will 
only consider trees of even sizes 2N where in the corresponding pure case, the reflexion vanish exactly for any flnite 
tree). It is instructive to write now the corresponding wave function as a function of the generation n, 

'4^pure{ 4) = 1 

'4^purei *^) ~ ^ 

'4^pure{ 2) — 1 

'4^pure{ 1) — ^ 

IppureiO) = 1 



V'pure (f ) 
'4'pure (2) 
'4^pure (3) 



^ 

'k 

1 

i 

1 



V'p«re(4) ... (23) 

This exponential decay of the wavcfunction in this dclocalizcd case is very peculiar to the tree geometry : it is imposed 
by the energy conservation and by the exponential growth of the number of sites with the generation n. 



C. Statistics over the disordered samples of the transmission Tjv 

1. Numerical pool method 

If one wishes to study numerically real trees, one is limited to rather small number of generations of order N^ax ~ 
12, 14 (sec for instance the study [37] based on exact diagonalization) because the number of sites and thus of random 
energies grows exponentially in N. From the point of view of convergence towards stable probability distributions via 
recursion relations, it is thus better to use the so-called 'pool method' that will allow us to study much larger number 
of generations (see below Eq. 27). The idea of the pool method is the following : at each generation, one keeps 
the same number Mpooi of random variables to represent probability distributions. Within our present framework, 
the probability distribution P„(-R) of the Riccati variables R at generation n will be represented by a pool of Mpooi 

complex values {Rn \ ■■, i?!*^''""'''}. It is convenient from now on to change the notation n 2N — n with respect to 
Fig. 1 so that n = now corresponds to the contacts with the outgoing wires (see Fig. 1). The numerical results 
presented below have been obtained with the following procedure : 

(i) the random variables of the initial pool are given by (Eq. 13) 

R<y{j) = -eo{j)-i (24) 

where e(j) are independent random variables drawn with Eq. 4. 

(ii) From the random variables Rn-iij) of the pool at generation n — 1, the pool at generation n is constructed as 
follows. For each j = 1, 2, ..Mpooi, one generates a new random energy e„(j) with the law of Eq. 4 and one draws K 
random indices {ii(j), ..jif (j)} among the pool of generation (n — 1) to construct the variable Rn{j) as follows(Eq. 
14) 

'^^-'""'-S/CT)) '''' 
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(iii) From these pools of Riccati variables, one may compute via Eqs 16 and 17 a pool of total transmission T„(j) 
for trees of n generations using 

T„(j) = l-|r„(i)|2 = l- 

For the Anderson model on the Bethe lattice, the pool method has been already used, in particular in [9] with a 
pool Mpooi = 1800 with a number N^ax ~ 30 of generations, in [17] with pools up to Mpooi = 10000, and in [20] with 
pools up to Mpooi = 16384 with a number N^ax ~ 400 of generations. The pool method is also very much used for 
disordered systems on hierarchical lattices (see for instance [38-40]), 

In the remaining of this section, we present the numerical results obtained with a pool of size 

Mpooi = 10^ with a number of generations N < N^ax = 34.10^ (27) 

We have also results for a pool of size Mpooi = 10^ with a number of generations N < Nmax = 24.10^ to see how the 

results change with the pool size. However, the number of generations Ng for this bigger pool Mpooi = 10^ has turned 
out to be less precise in the critical region. All figures shown below thus corresponds to data obtained with the pool 
of size Mpooi = 10^. 

As is usual with the pool method [40] , the location of the critical point depends on the pool, i.e. on the discrete 
sampling with Mpooi values of probability distributions. It is expected to converge towards the thermodynamic critical 
point only in the limit Mpooi oo (see the discussion of section B3 in Appendix B). Nevertheless, for each given 
pool, the critical behaviors with respect to this pool-dependent critical point usually allows a good measure of critical 
exponents [40]. For instance, for the pool of size Mpooi = 10^, the critical value Wc of the disorder strength where 
the localization-delocalization occurs for the Landauer transmission is of order 

WciMpooi = 10^) ~ 16.99.. (28) 

For the pool of size Mpooi = 10^, we find that it is higher and of order Wc{Mpooi = 10^') ~ 17.32... This rather 
important shift of the pseudo-critical point with the pool size which has already been seen in [17], seems to be due to 
the very slow logarithmic convergence of traveling wave velocity in the presence of cut-off (see the discussion of section 
B 3 in Appendix B). So we stress that here, in contrast to [9, 17], our goal is not to determine the true thermodynamic 
mobility edge Wc{+oo), but instead to understand the critical behaviors of the finite-pool results with respect to the 
pool-dependent critical point of Eq. 28. 

We will first discuss the behavior of the typical transmission Tj^^ defined by 

ln(T^^^') = h^ (29) 

as a function of the number A'' of generations and disorder strength W, before we turn to the distribution around this 
typical value. 



i + RnU) 



RnU) 



(26) 



2. Exponential decay of the typical transmission Tj^^ in the localized phase W > Wc 

In the localized phase, one expects that the typical transmission defined by Eq. 29 decays exponentially with 
the number N of generations 

where ^loc represents the localization length that diverges at the delocalization transition 

^ioc{W) {W-Wc)-"'"^ (31) 

We show on Fig. 2 our numerical results for the pool of size Mpooi = 10^ (Eq. 27) : the exponential decay with 
N of Eq. 30 is shown on Fig. 2 (a) for various disorder strength TV. The corresponding slope l/S,ioc{W) is shown as 
a function of W on Fig. 2 (b) : we find that this slope vanishes linearly in {W — Wc), in agreement with the exact 
result [10, 11] 

i^ioc = 1 (32) 
and in agreement with Figure 7 of Ref. [20] concerning a random-scattering model on the Cayley tree. 
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FIG. 2: Exponential decay of the typical transmission T^'' = e'""^" in the localized phase W > Wc ■ (a) Linear decay of InTjv 
as a function of the number A'^ of generations (see Eq. 30). (b) Behavior of the slope l/£,ioc(W) as a function of the disorder 
strength W : it vanishes linearly l/^ioc{W) ~ (V - wr°T""' with WT"' ~ 16.99 (see Eq. 28 ) and uioc = 1 (see Eq. 31) 



3. Behavior of the typical transmission in the delocalized phase W < Wc near criticality 





FIG. 3: Behavior of the typical transmission T^'' of the infinite tree in the delocalized phase : (a) InT^'' = \nToo{W < Wc) 
as a function of the disorder strength W (b) same data in a log-log plot to measure the exponent of the essential singularity of 
Eq. 34 : ln(- InT^^P) as a function of ln{Wc - W) : the slope is of order kt 0.25 



In the delocalized phase, the typical transmission remains finite in the limit where the number of generations N 
diverges 



\nTN{W < Wc,N) 



N~ 



In Too (W^ < Wc) > -oo 



As shown on Fig. 3, we measure the following essential singularity behavior of the typical transmission 

\nT^(W < Wc) ~ -{Wc-W)-''^ with 



W~*Wa 



(33) 



(34) 



The presence of essential singularities in transport properties near the localization transition on the Bethe lattice 
has been found in [11] via the supersymmetric formalism (see also [12-15] where similar results have been obtained 
for the case where Anderson tight-binding model is replaced by a non-linear a-model). In particular, Eq. 71 of Ref. 

[11]) states that the leading critical behavior of the diffusion constant is given by : InD \E - Ec\-^^^ (Eq. 71 of 

Ref. [11]). We note that the exponent in this essential singularity is 1/2 instead of the exponent of order 1/4 that 
we measure. The reason for this difference could be that the results of Ref. [11] are based on the computation of the 



8 



disorder-averaged two-point density- density correlation function (sec Eq. 2 of Ref. [11]), whereas our numerical results 
concern the typical value of Landauer transmission, and not the disorder- averaged transmission which is expected 
to be governed by rare events (see below the section II C 5). Also it is not clear to us what exactly represents the 
diffusion constant computed in [11], because it is known that on the Bethe lattice the dynamics is not diffusive but 
ballistic in the delocalized phase [41], and that more generally the random walk on the Bethe lattice is not diffusive 
because the tree-geometry induces an effective bias away from the origin (see for instance [42] and references therein). 
As a consequence, some implicit reinterpretation of the Bethe lattice seems to underlie the statements of Ref. [11] 
that makes difficult a precise comparison with our present numerical results. 



4- Finite-size scaling in the critical region 



If there exists some finite-size scaling in the critical region for the typical Landauer transmission of the form 

lnTjv(W") -NP-^G (^N^/^^^iWc - W)^ (35) 

the matching of our results in the localized phase (see Eq. 30) and in the delocalized phase (Eq. 34) requires a 
finite-size correlation length exponent v^^ of order 

l^T^ = noc + KT = l + KT^ 1.25 (36) 

In another traveling/non-traveling phase transition studied in Ref [30] (see the summary in Appendix C), it has been 
obtained that the finite-size scaling exponent vps is determined by the relaxation rate towards the finite value in the 
non-traveling phase. We have thus studied the relaxation length towards the finite value in the delocalized phase. We 
find that our data for In Tjv are compatible with the form 



lnTjv(VF < Wc) =i -{Wc - Wy-^ON (37) 
where ajv is a random stationary process as a function of N . We find that its autocorrelation function is exponential 

C{N) ~ e"«re!a"(W) (38) 

and we measure that the relaxation length ^reiax{W) diverges with an exponent 

^relax{W)cX-^^^—^^^y^ with Urelax ^ 1-21 (39) 

of the order of the exponent u^^ of Eq. 36. We thus obtain that the critical properties are qualitatively similar to the 
critical properties described in Ref [30] (see the summary in Appendix C) : the traveling phase is characterized by a 
velocity that vanishes linearly, but the finite size scaling is governed by the relaxation length towards the asymptotic 
finite value of the non-traveling phase. Exactly at criticality, we thus expects the following stretched exponential 
decay of the typical transmission 



InTjv(Wc) -NP-^ (40) 

where the exponent pT is related to the other exponents by (see the scaling relations of Eqs C9 and Cll in Appendix 
C) 

= ^ = 1 - -is (41) 
Prom our previous estimate of the exponent kt — 0.25, this would correspond to the numerical value 

PT = ^ 0.2 (42) 



We have not been able to measure this stretched exponential behavior exactly at criticality from our data, because a 
precise measure of the exponent pt would require to be exactly at the critical point. 
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InT. 



FIG. 4: Evolution with N of the probability distribution Pjv(lnTjv) of the logarithm of Landauer transmission : (a) In 
the delocalized phase (here W = 10), the probability distribution Pjv(lnTjv) does not move with A^, and vanishes with a 
discontinuity at the boundary InTjv — 0. (b) In the localized phase (here W — 25), the probability distribution Pjv(lnrjv) 
moves with A'^ as a traveling wave of fixed shape. 



5. Distribution of the logarithm of the transmission 

Up to now, we have only discussed the behavior of the typical transmission of Eq. 29 as a function of A'' and W. 
We now turn to the probability distribution of the logarithm of the transmission around its averaged value, i.e. we 
consider the distribution of the relative variable 



= InT, 



N 



InT, 



N 



(43) 



We find that as oo, this variable remains finite not only in the delocalized phase where In Too is finite, but also in 
the localized phase where In Tjv decays linearly in N (Eq. 30). This means that in the localized phase, the probability 
distributions P;v(liiT/v) actually propagates as a traveling wave with a fixed shape around its moving center InT/v as 
shown on Fig. 4 (b) . This phenomenon has already been seen for random scattering models on the Bethe lattice (see 
Fig. 8b of Ref. [20]). This is in contrast to the broadening with L observed in low dimensions d ~ 1,2,3 (see Eq. 2) 
: the Cayley tree thus corresponds to oJcayiey = in Eq. 2. 
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FIG. 5: (a) Logarithm of the probability distribution P{u) of u = In Tat — InTjv as W varies in the localized phase for various 
disorder strengths W = 17,25,50,100 : the exponential decay as e"''" corresponds to slopes of order I3(W = 17) ~ 0.48, 
I3{W = 25) ~ 0. 43, I3{W = 50) ~ 0.35, I3{W = 100) ~ 0.33. (b) Width a = (y f^'^ of the probability distribution of 
u = In Tn — In Tjv as a function oiW: it remains finite both in the localized phase and in the delocalized phase, but it presents 
a cusp singularity at criticality. 



The fact that the shape is fixed can be used numerically to measure more precisely the tails of this probability 
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distribution by accumulating data over iterations : we show on Fig. 5 (a) the histograms of u obtained by this 
procedure for various disorder strength W. An essential property of this distribution is the right exponential tail 

Pw{u) e-'^(^)" (44) 

In terms of the rescaled transmission t = Tn/T^^ = (see Eq. 43), this corresponds to the power-law decay 

P(i = TWTr),J-^^ (45) 

We find that the selected exponent f3seiec{W) slightly grows as the disorder strength W decreases, from a value of 
order P{W = 100) ~ 0.33 for the strong disorder W ~ 100 towards a value of order P{W = 17) ~ 0.48 near criticality 
(see Figure 5 a). The velocity of the traveling wave propagation of the whole distribution is actually determined by 
this power-law tail, as explained in detail in Appendix A. In particular, the critical value (^{W Wc) = 1/2 has 
been predicted in [9] (sec Eq. A33 in Appendix A) . An important consequence of the power-law tail of Eq. 45 with 
(3{W) < 1/2 is that all integer moments of the Landauer transmission will be governed by rare events. 

Finally, we show on Fig. 5 (b) the width of the probability distribution of In T/v as a function of : it remains 
finite both in the localized phase and in the delocalized phase, but it presents a cusp singularity at criticality because 
in the delocalized phase, the distribution presents a singularity at finite distance from the typical value, as is clearly 
visible on Fig. 4 (a) : the transmission is bounded by T < 1, i.e. the histogram of InTjv presents a discontinuity at 
InTiv = 0. 



III. STATISTICS OF EIGENSTATES 



The Landauer transmission studied in the previous section is of course the most appropriate observable to charac- 
terize the transport properties and to find the transition between the conducting/non-conducting phases. However, 
one expects that these transmission properties that emerge when the disordered sample is linked to incoming and 
outgoing wires are related to the nature of eigenstates of the disordered sample in the absence of these external wires 
(see Fig. 6) To determine whether a normalized eigenstate 

is localized or delocalized, the usual parameters are the Inverse Participation Ratios (I.P.R.) 

X 

of arbitrary power q. Another important quantity to characterize the spatial extent of an eigenstate is its entropy 

S^-Y,\i,{x)\Hn\i,{x)\^ (48) 

X 

The normalization condition yields the idcndity Iq^i = 1, and the entropy corresponds to 

5 = -ag/g|g=l =-5gln/g|g=l (49) 

On an hypercubic lattice of size L in dimension d, containing Ml = L'^ sites, the delocalized phase is characterized 
by the following typical behaviors 

hiTg ~ -d(g- l)lnL = -(g- l)lnA/'L (50) 

and 

S ^d\nL = \nML (51) 

These scalings arc the same for pure homogeneous eigenfunctions with |V-'(2;)| = l/£'' for all sites x G L'*. At criticality, 
the scalings of the I.P.R. involve a whole series of non-trivial exponents and the wave function is said to be multifractal 
(see the review [7]). Finally in the localized phase, the I.P.R. and the entropy are finite. 

In the present section, we discuss the statistical properties of these I.P.R. and entropy for zero-energy eigenstates 
for the Cayley tree of Fig. 6 as a function of the number N of tree generations and of the disorder strength W . 
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FIG. 6: Caylcy tree with branching ratio K = 2 where each interior site has /sT + 1 = 3 neighbors. On the figure the tree ends 
at generation 2N = 4. 



A. Reminder on the Miller-Derrida framework to construct eigenstates 

We refer to [17] where it is explained how eigenstates of finite trees can be constructed and how the density of states 
can be obtained. Note that on a tree, boundary sites dominate so that one needs a substraction procedure to obtain 
the appropriate bulk density of states. In the following, we are not interested into the density of states, which does 
not contain any information on the localized/delocalized nature of the spectrum. Wc wish to study instead the spatial 
properties of eigenstates of zero-energy E = (center of the band). The Schrodingcr equation 9 yields as before the 
recursion of Eq. 14 for the Riccati variables defined in Eq. 11. The difference with the scattering case is now in the 
boundary conditions 

R{2N,j) = -e{2N,j) (52) 

that replace Eq. 13. As a consequence, the Riccati variables are now real (and not complex). The energy E ^ will 
indeed be an eigenstate only if the Schrodinger equation is also satisfied at the center of the finite tree at the center 
site 

K+l 

Since the on-site energy e(0) is a random variable drawn with some distribution, we may consider that we choose e(0) 
to satisfy Eq. 53 to obtain a typical eigenstate of zero energy and to study its spatial properties. Let us first describe 
zero-energy eigenstate in the pure case to stress the peculiarities of the tree geometry. 

B. Structure of zero-energy eigenstate on the pure tree 
In the pure case where all on-site energies vanish e{n,j) = 0, the zero-energy eigenstate with radial symmetry reads 

V'2n-l = (54) 
V'2n =V'o(-^) (55) 



where the root amplitude ipo is determined by the normalization condition 

1 = J2 IV'WP = IV'oP + iK+ l)i^|V2|' + (K + l)i^3|V'4p + ... + {K + l)ir2iv- V2iV 



IV'ol 



(56) 
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This state is the analog of the scattering state described in section TIB. In particular, one sees again an exponential 
decay of the wavefunction in this delocalized case which is very peculiar to the tree geometry : in the normalization 
condition of Eq. 56, this exponential decay is exactly compensated by the exponential growth of the number of sites 
at generation 2n, so that all generations 2n = 2,4, ..2N carry exactly the same weight in the normalization condition, 
and the root weight I'f/'oP vanishes as 1/A'' in the thermodynamic limit, where (2A'') is the number of generations of 
the tree. 

For this state, the usual Inverse Participation Ratio of Eq. 47 for g = 2 reads 

ir'{2N) = lV'(a;)l' = IV'ol' + {K+ l)K\ij^\^ + {K+ l)K^\M'' + - + {K + 1)1^^^- V2iv|' (57) 

K+l l-(i.) 



1^0 1' 



1 + 



1 - A 



(58) 



More generally for g > 1, one obtains the following decay with the number of generations 



In terms of the total number of sites 

2N-1 



l^ol^'-^ for g>l (59) 



N{2N) =l + (if + l) ^ i^" = l + (i^ + l)^^__ = L^^-J (60) 

this corresponds to 



K-l K 



ln7P"''(2A/-) oc -q\nN cx -q\n{\nN{2N)) for q>l (61) 

N—>-oo N—^oo 

This behavior is thus very anomalous with respect to the corresponding decay as {—{q — l)lnjVi,) for the I. P.P. of 
pure states in finite dimension d (see Eq. 50). 

Moreover, since the eigenfunction normalization yields the idendity 7g=i = 1, one sees that there exists some 
discontinuity as g ^ 1. In particular, the entropy of Eq. 48 which can usually be computed as the derivative of Eq. 
49 should be computed directly here, and one obtains 

5P«r(27V) = -|V(0)|2 In \^{0)f - {K + l)K\i;{2)f In \ip{2)\^ - ... - {K + l)K^^-^\iP{2N)f In \i}{2N)f (62) 



1 + ^^4^n) In |^(0)|2 + ^^^N{N + l)lnK 
K J K 



(63) 



= 1^(0)1 

Using Eq 56, one obtains the following behavior for large N 

( K + \ \ i^iV(7V+ l)lnX 
5P«r = In 1 + + ^ , ^ r,^' . ^ NlnK 64 

In terms of the total number of sites of Eq. 60, the entropy is proportional to the logarithm of the number of sites 

SP^r^2N) ~ NlnK ~ -lnAf{2N) (65) 

iV— »oo N^oo 2 

that should be compared with the growth as lnM{2N) for pure eigenstates in dimension d (see Eq. 51). 

In conclusion, the tree geometry induces very anomalous scalings in (In In A/") for the logarithm of I.P.R. of pure 
states (see Eq. 61) with respect to the case of finite dimension d, whereas the entropy is well behaved in (InA/") (see 
Eq. 65). In the disordered case, we expect to observe similar behaviors, whereas the I.P.R. and the entropy will 
remain finite in the localized phase. 

C. Recursion relation for Inverse Participation Ratios (I.P.R.) 

To compute the Inverse Participation Ratio Iq (Eq. 47) via recursion, one needs to introduce besides the Riccati 
variable of Eq. 11 the auxiliary variables defined by 

ip(anc[n,j)) ^ ^p[anc[n,j)j 
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where the sum is over ah sites / that are descendants of the site {n,j) of the tree. These variables satisfy the following 
recurrence inside the tree 1 < n < 2N — 1 



1 



i?2(n,j) 



K 



and the initial conditions at the boundaries 

C(«)(27V,j) = 

At the central root, on needs to impose the normalization 

i = ^|V'(x)p = |^(o)p 



m=l 



R^{2N,j) 



K+l 



l+^CW(l,m) 



m=l 



(67) 



(68) 



(69) 



that determines the weight |V'(0)P of the root in terms of the variables C^^^ of the branches. 

To compute the inverse participation ratio of parameter q of Eq. 47, one needs the variables C^'^ together with the 
variables C'^' 



/, = ^|V^(a;)|2« = 1^(0)1 



2q 



K+l 



1+ J2 c^'Hi^H 



K+l 

l+^C(«)(l,m) 

m—l 



K+l 



1 1 



(70) 



1+ ^ C^^\l,m) 



D. Recursion for the entropy 



To compute recursively the entropy of Eq. 48, one needs similarly to introduce the auxiliary variable 



where the sum over I denotes the sum over the site (n,j) and all its descendants : a{n,j) thus represents the entropy 
for the branch containing (n, j) and its descendants. 
The initial conditions at the boundaries are simply 

a{2N,j)=0 (72) 

and the recursion inside the tree 1 < n < 2N — 1 can be written as 

/ .^ Em=i <^^^^ {desrnjn, j))cr{deSm{n, j)) 

where the first term represents the weighted contribution of the branches entropies and where the second term 
represents the mixing entropy 

K 



with the weights 



-Polnpo- X! Pni''iT^Pm' (74) 



Po = f? (75) 



i^ + El=iC('KdeSm{n,j))) 

(77) 



normalized to po + Ylm'=iPm' = 1- At the central site of the tree, one uses the same formula but with {K + 1) 
branches instead of K branches. 
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E. Numerical results on the statistics of eigenstates 



1. Numerical pool method 



To study the statistical properties of eigenstates, we have used again the pool method explained in II CI. The 
only difference is that for the transmission, we have followed the recursions for the complex Riccati variables, whereas 
here we follow the recursions for the real Riccati variable R and for the auxiliary variables C'^-*, C'^^-' and a described 
above. 



2. Statistics of the entropy of an eigenstate 




FIG. 7: Statistics of the entropy S{N) = — In |i/)(a;)|^ of a normalized zero-energy eigenstate : (a) Evolution of the 

typical entropy Styp{N) = e''^'^'^' with the number A'^ of generations : it grows linearly in A'^ in the delocalized phase, whereas 
it remains finite in the localized phase, (b) Asymptotic width a — (u^)^/^ of the relative variable u = In S — In S as a function 
of the disorder strength W : the width converges to (as 1 /TV) in the delocalized phase, whereas it is finite in the localized 
phase 



We first consider how the eigenstate entropy defined in Eq. 48 evolves with the number N of generations. As shown 
on Fig. 7 (a), the typical value grows linearly in N in the delocalized phase 



Styp{N, W <W,) = e 



\nS(N) 



aN 



(78) 



where the factor a varies smoothly with W and does not vanish continuously near the critical point. In this delocalized 
phase, we moreover find that the width of the relative variable v = In 5jv — In Sn decays to zero as N oo 



(In Sm — In SnY'' 



1/2 1 
W<W^ N~Qo N 



(79) 



In the localized phase on the contrary, the typical value and the width of the relative variable v = In Sn — In Sn 
remain finite as ^ +oo (see Fig. 7 (a) and (b)). 

As shown on Fig. 8, we find that the typical entropy of an infinite tree of the localized phase diverges with the 
following power law near criticality 

The entropy measures the size of the region where the weight |-0(a:)|^ is concentrated, whereas the localization length 
Cioc ~ ^/{W — WcY'""^^ describes the far exponential decay of the transmission T (Eq. 30). We thus conclude that 
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in the localized phase, besides the localization length ^/oc ~ ^/{W — Wc)'''°=~^ known since [10], there exists a larger 
diverging length scale 

('(^>^'\„-„„w^m^' (81) 

that characterizes the size where the weight | ■0(2^)1^ is concentrated. 



3. Statistics of the root weight \tpN{0)f and of the I.P.R. I2 

As explained above, the root weight \iPn{0)\'^ and of the I.P.R. I2 are determined by the auxiliary variables C*^*^ 
introduced in Eq. 66. We thus expect that the localized and delocalized phases correspond to the following behaviors 
for the auxiliary variables C''^ 

(i) in the localized phase, the auxiliary variables C*^''-' will remain finite random variables as the number of generations 
diverge TV — > (X). Then the root weight |^Ar(0)p and the I.P.R. Ig remain finite as ^ cxo. 

(ii) in the delocalized phase, the auxiliary variables C^'^^ will instead grow exponentially with the number N of 
generations with some Lyapunov exponents Ag > defined by 

lnC(9) ~ XgN (82) 

The corresponding typical behaviors of the root weight and of the I.P.R. then reads 

ln|Vjv(0)|2 ~ -AiTV (83) 

iV— >(x) 

and 



In/ ~ _(gAi-A,)iV (84) 

We find numerically that the probability distributions of auxiliary variables (InC^) move as traveling waves of 
velocity Xq (see Eq. 82) with fixed shape (see Figure 12 of Appendix B where the fixed shape is shown for g = 1 and 
q ~ 2). We explain in Appendix B how the Lyapunov exponents A, can be determined via a tail analysis that yields 
the identity (see Eq. B13) 

Xq = qXi (85) 

so that the logarithm of the I.P.R. Iq decays slower than linearly in N , because the coefficient in Eq. 84 exactly 
vanishes. Physically, one could expect the decay to be of order (InA^) as in Eq. 61 concerning pure states on the 
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Cayley tree, but we have not been able to measure the behavior of the I.P.R., because within the pool method that 
we use, it turns out that the pool-dependent critical point is not the same for the variables C^'^-' and for C*-^-* (see 
more details in section B 3), so that it does not seem easy with the pool method to extract reliable results concerning 
h- 

In the remaining of this section, we thus focus on the statistical properties of the root weight |-!/'(0)p that is 
determined by the auxiliary variable C*^^^ via Eq. 69. The fact that the pool-dependent critical point for C^^^ is 
the same as the critical point found for the Landauer transmission can be understood from the tail analysis of the 
Appendices A and B that involve exactly the same integral kernel. 




FIG. 9: (a) In the delocalized phase, the logarithm of the root weight decays linearly ln|i/;jv(0)p ~ —\\{W)N : the figure 
shows that the slope \i{W) vanishes linearly \i{W) ~ {Wr°^ - W) (see Eq. 86 with ui^^ = 1). (b) Width a = (z^f^^ of the 
relative variable z = \n |i/)]v(0)p — In |?/)jv(0)P in the limit A*' -f oo as a function of the disorder strength W : it remains finite 
both in the localized and delocalized phase, but it presents a cusp at the critical point Wc- 



In the delocalized phase, the root weight |7/;Ar(0)p presents the typical decay of Eq. 83. We show on Fig. 9 (a) that 
the corresponding Lyapunov exponent Xi(W) vanishes linearly in the critical region 

Xi{W) ~ (Wc-W) (86) 

i.e. we find the same exponent vioc = 1 as for the divergence of the correlation length of the localized phase (see Eq. 
32). This can be understood from the tail analysis presented in the Appendices A and B. 

We find that the width of the relative variable z = In |'07v(O)P — In |-!/)Ar(0)p remains finite in the limit N +oo 
both in the localized phase and in the delocalized phase : its behavior as a function of the disorder strength W 
presents a cusp at Wc as shown on Fig. 9 (b). 




FIG. 10: Critical behavior of the root weight |i/)oo(0)p of an infinite tree in the localized phase : (a) In |'!/'cx)(0)P as a function of 
the disorder strength W (b) ln(— In IV'oo(O)p) as a function of ln(W — Wc) to measure the exponent of the essential singularity 
of Eq. 87 : we measure asymptotically a slope of order kq ~ 1- 
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In the localized phase, the root weight remains finite as -/V — > oo. As show on Fig. 10, we measure the following 
essential singularity 

1 



^"^-^'^ wZv, - {W -W.ro -tl^ -0 - 1 (87) 

As in section II C 4, wc now discuss the finite-size scaling in the critical region. If there exists some finite-size scaling 
in the critical region for the typical root weight, the matching of our results in the dclocalized phase (see Eq. 83 and 
86) and in the localized phase (Eq. 87) requires a finite-size correlation length exponent i'q^ of order 



,FS 



1 + Ko ^ 2. (88) 



We have performed the analysis described in section IIC4 and we measure that the relaxation length ^reiax{W) 
towards the finite value of Eq. 87 diverges with an exponent 

Pretax (W) OC j^^^^—^^y-^ with l/relax - (89) 

of the order of the exponent Vq^ of Eq. 88. We thus obtain that the critical properties are again qualitatively 
similar to the critical properties described in Ref [30] (see the summary in Appendix C) : the traveling phase is 
characterized by a velocity that vanishes linearly, but the finite size scaling is governed by the relaxation length 
towards the asymptotic finite value of the non-traveling phase. Exactly at criticality, we thus expects the following 
stretched exponential decay of the typical transmission 



ln|V'jv(0)|2|iy=iy^ --Af"" (90) 

where the exponent po is related to the other exponents by (see the scaling relations of Eqs C9 and Cll in Appendix 
C) 

^0 = -FS = 1 - ^ (91) 
Prom our previous estimate of the exponent kq — 1-, this corresponds to the numerical value 

(92) 

Again, as explained after Eq. 42, we have not been able to measure this stretched exponential bc^havior exactly at 
criticality from our data, because a precise measure of the exponent po would require to be exactly at the critical 
point. 



IV. CONCLUSIONS AND PERSPECTIVES 



In summary, for the Anderson model on the Bethe lattice, we have studied numerically the statistics of the Landauer 
transmission Tjv and the statistics of cigcnstatcs at the center of the band i? = 0, as a function of the disorder strength 
W and the number N of generations. We have shown that both the localized phase and the delocalized phase are 
characterized by the traveling wave propagation of some probability distributions. In the text, we have presented 
detailed numerical results, and in Appendices A and B, wc have explained how the velocities of the traveling waves 
are determined by the tails, via the properties of the integral kernel introduced in [9] . The Anderson transition then 
corresponds to a traveling/non- traveling critical point for these traveling waves, and the critical properties obtained 
are very similar to the traveling- wave phase transition studied in [30] : (i) the finite value of the non-traveling phase 
presents an essential singularity (ii) the relaxation length towards this essential singularity determines the finite-size 
scaling in the critical region (iii) the finite-size correlation length exponent i'fs is different from the value z/joc = 1 that 
govern the vanishing of the velocity in the traveling-wave phase. We thus hope that in the future, these properties 
for Anderson localization on the Cayley tree will be better understood by extending the methods that have been 
developed recently for the class of the Fishcr-KPP traveling waves (see [30, 43] and references therein). 

How results obtained on the Bethe lattice are related to properties of Anderson localization in finite dimension d 
is of course a difficult question. During the history of localization, many values for the upper critical dimension have 
been proposed such as = 4,6,8, -l-oo. Even within the supersymmetric community, there seems to be different 
interpretations. Mirlin and Fyodorov [45] consider that the essential singularities that appear on the Bethe lattice. 
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arc directly related to the exponential growth of the number of sites with the distance, and will become conventional 
power-law behaviors as soon as d < cxd, so that the upper critical dimension is actually = +00. On the contrary, 
Efetov [15] argues that the results obtained for the Bethe lattice are relevant to Anderson transition in high dimension 
d when reinterpreted within the so-called 'effective medium approximation'. For the Landauer transmission, the upper 
critical dimension dc can be defined as the dimension where the exponent u{d) concerning the width of the sample- to- 
sample distribution in the localized phase (see Eq.2) vanishes ui{dc) = 0. This means that for d > dc, the probability 
distribution would travel as a traveling wave in the whole localized phase. To determine whether dc = +00 or 
dc < +00, it would be thus interesting to understand whether the traveling wave propagation with a fixed shape is 
possible only on trees or whether it can also occur in sufficiently high dimension d. 

In low dimensions, the probability distribution of the logarithm of the transmission is known to broaden with 
some exponent u{d) > (see Eq. 2) in the localized phase, whereas the probability distribution of the logarithm 
of I. P.P. is known to shrink with the system size in the delocalized phase [46] (see also Fig. (4a) of [49] where the 
same phenomenon occurs for the directed polymer in 1 + 3 dimensions). However traveling wave propagations of the 
probability distribution P{\nlq) of Inverse Participation Ratios Iq have actually been found in finite dimensions but 
only exactly at criticality. This phenomenon has been first obtained for the power-law random band matrix mode l [47] 
and has been then observed for the Anderson model in dimensions d = 3, 4 [48] : the motions of the typical values In/^ 
with (InL) determine the multifractal spectrum, whereas the reduced variables y = Iq/Iq^^ keep fixed distributions 
presenting power-law tails with g-dependent exponents. We have actually observed the same behavior for the directed 
polymer in 1 -|- 3 dimensions exactly at the localization/delocalization transition (see Fig. (3a) and (3b) of Ref. [49]). 
In all these cases, it would be thus very nice to better understand the relations between the tails of the distributions 
and the motions of the typical values, and to identify the wave equation that underlies these traveling waves. 

Finally, our work raises once again the question whether Anderson transitions are characterized by several length 
scales that diverge with various v exponents. Here for the Cayley tree, besides the exponent i>ioc = 1 that has 
been predicted for a long time [10], we have found other diverging lengths with the following physical meanings, 
(i) in the localized phase, besides the exponent vioc = 1 associated to the divergence of localization length S^ioc ~ 
1/{W — WcY^""^^ that describes the far exponential decay of the transmission T (Eq. 30), we have found that the 
critical behavior of the entropy is governed by another diverging length scale with z/g ~ 1.5 (Eq. 81) that characterizes 
the size where the weight |?,i'(.T)p is concentrated, (ii) in the delocalized phase, we have found that the Landauer 
transmission reaches its asymptotic value after a length diverging as v^^ ~ 1.25 (Eq. 36) that governs the finite-size 
scaling of the transmission in the critical region. Even if the precision of numerically measured critical exponents can 
be always discussed, we feel nevertheless that our numerical results arc not compatible with the existence of the single 
exponent uioc = 1- Note that for the directed polymer on the Cayley tree, we have also found previously that the 
critical properties involve two exponents v = 2 and v' = 1 [50]. Whether there exists various exponents v for Anderson 
model in finite dimension d has been debated for a long time. The majority of papers on Anderson localization seem 
in favor of a single v (see the reviews [2-7]), but it seems to us that numerical papers actually study always the 
same type of observables. In particular, we have not been able to find numerical studies concerning the entropy of 
eigenstates. On the theoretical side, we are not aware of many theories in favor of various v exponents, except (i) 
the supersymmetric studies of Refs [14, 15] that involve two different diverging length scales, called respectively the 
localization length and the phase coherence length, (ii) the pseudo-delocalization transition of the random hopping 
model in dimension d = 1, that is characterized by the two exactly known correlation lengths exponents v = 1 and 
V = 2 (see [51] and references therein), as in the other models described by the same strong disorder fixed point (see the 
review [52]). We feel that if localization models exhibit several v exponents both on the Cayley tree which represents 
some mean-field d = qg limit and in some models in dimension d = 1, the possibility of various v for Anderson 
transition in dimension d = 3 should be reconsidered by studying in detail the statistics of various observables. 

APPENDIX A: TAIL ANALYSIS FOR THE TRAVELING WAVE OF THE LANDAUER 
TRANSMISSION IN THE LOCALIZED PHASE 

In this Appendix, we translate the analysis of Ref. [9] concerning the distribution of the self-energy in the localized 
phase for the traveling-wave propagation of the Landauer transmission discussed in section II. 

The recursion of Eq. 14 for the complex Riccati variable Rn reads more explicitly in terms of its real and imaginary 
parts Rn = Xn — iYn with X„ e] — 00, -|-oo[ and Yn € [0, -|-oo[. 
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Yn 



K 



Xn-i{m) 



K 

E 



Yn-i{m) 



'^Xl_,{m) + Y^_,{m) 



(Al) 
(A2) 



where the random on-site energies e„ are drawn from the flat distribution of Eq. 4. The Landauer transmission of 
Eq. 26 then reads 



T„ = 1 - 



x„ + i(i-r„) 



-Xn + i{l + Yn) 



4y„ 



XI + {Yn + 1)2 



(A3) 



In the dclocaUzcd phase, this recursion of Eq. A2 is difficult to analyze because one has to find the joint distribution 
of the two finite variables [X, Y) that remain stable upon iteration. In the localized phase however, the problem is 
simpler [9] as we now recall. 



1. Linearized recursions in the localized phcise 

In the localized phase, if the imaginary part Yn converge towards zero exponentially in n as 

^ n — ^ yn 



(A4) 



where v > and where y„ remains a finite random variable upon iteration, one has then to study the simpler 
recurrence for large n [9, 17] 



Xn 



e Vn 



y ' 



K 



= E 



y„_i(m) 



Note that the form of Eq. A4 corresponds to a traveling wave of velocity v for the variable (Inl^) 

lny„ = -un + lnt/n 



(A5) 



(A6) 



The real part X„ satisfies now a closed recurrence independent of the Yn-, and its stable distribution P*(X) satisfies 
the closed equation 



P*{X) 



I dep{e) J < 



dXiP*(Xi)... / dXmP*{X„,)S 



X + e- 



^ 1 



m=l 



An important property of this distribution is that it presents the power-law tail 

KP*{0) 



P*{X) 



X^±oo X"^ 



(A7) 



(A8) 



because whenever one of the K variables on the right handside of Eq. AT is close to 0, the variable X of the 

left-handsidc is large with X ^ —1/Xi. 

The stable joint distribution P* {X, y) satisfies 



P*{X,y)^ Jdepie) j dX^dy^P*{Xuy^)... j dXn,dymP*{Xm,ym)5 ^ + e+E^ 



K 

E 

m=l 



e-^Xl 



(A9) 
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Tail analysis 



The idea of [9] is to look for the power-law tail in y that is compatible with the recursion equation, i.e. one assumes 

A{X) 



P*{X,y) 



J/— »(x) y 



1+/3 



(AlO) 



with some exponent < /3 < 1. Note that, in the traveling wave language of Eq. A6, this is equivalent to look for 
solutions with an exponential tail e"''" for the variable u = lny„ + nv = lny„. So this corresponds to the usual 
exponential tail analysis of fronts [28, 29], except for the following difference : in usual studies of propagation into 
unstable phases, it is the 'forward tail' that has an exponential decay and that determines the velocity, whereas in 
our present case, it is the 'backward tail' that determines the propagation. 

In Laplace transform with respect to y, the power-law decay of Eq. AlO corresponds to the following singular 
expansion near the origin 



P + OC P + OO 

P*{X;s) = J^ dye-'yp*{X,y) =P*{X)-I dy{\ - e-'y)P*{X,y) ^ P*{X) 



r+°° dv 

Jo s 



P''{X)-s''A{X) I 



+~ , (l-e-") 



{l-e~^)P*{X,-) 



(All) 



Equation A9 becomes for the Laplace transform P*{X; s) 
P*{X;s)-- 



dep{e) I dX^P*{X,,^)...J dXmP*{X^,j^)6 



X + e- 



^ 1 



m=l 



Using the expansion of Eq. All, one obtains that the function A{X) has to satisfy the eigenvalue equation 

e-^^A{X) =kJ ^Q{X + ^)^(^i) 
where the function Q{u) represents the stationary distribution of the variable u = —e — X^^=2 X~ 



(A12) 



(A13) 



Q{u) = J dep{e) J dX2P*{X2)... j dX^P*{Xm)5 
(see Eqs 6.5 and 6.6 in [9]). 



^ 1 



m^2 



(A14) 



3. Eigenvalue problem for an integral kernel 

The tail analysis thus leads to the eigenvalue problem the integral kernel appearing in Eq. A13 

KA{X) = K j ^Q(X + ^)^(^i) (A15) 

For each /? and W, the integral kernel is positive, and thus one expects some continuous analog of the Perron- 
Proebenius theorem : the iteration of the integral kernel will converge towards a positive eigenvector Ao{X) that is 
associated to the maximal eigenvalue Ao(/3, W) of the kernel. For instance, for the special case /? = 0, the solution is 
simply A{X) = P*{X) (see Eqs A7 and A14) and 

Ao(/3 = 0, VK) = X (A16) 

Except for this case /? = 0, we are not aware of any explicit solution for the eigenvalue Ao(/3, W) (even for the simpler 
case where the distribution of the on-site energies p{e) is a Cauchy law and where the distributions P^X) and Q{u) 
are also Cauchy laws [9, 17]). With Eqs A13 and A15, one concludes that at fixed W, each mode /3 is associated to 
the velocity v((3, W) ( Eq. A4) 



v{f3,W) = -^lnAo{l3,W) 



(A17) 



21 



4. Selection of the tail exponent and of the velocity of the traveling wave 

The selection of the tail exponent /3 of Eq. AlO and of the corresponding velocity of Eq A4 usually depend on 
the form of the initial condition [26, 28, 29]. In our present case, the initial condition is completely localized (see Eq. 
13) 

P,miY) = S{Y - 1) (A18) 

In this case, one expects that the solution that wiU be dynamically selected [26, 28, 29] corresponds to the tail 
exponent PseieciW) and to the velocity VgeiedW) = v{f3seiec{W),W) determined by the following extremization 



4lnAo(/J,l^)-r^^°^''^^ 



(A19) 

/3=/3seiec("') 



^ p Ao{P,W) 
The critical point is then determined by the two conditions 

v{p,,W,) =0 (A20) 
di^XPcW,) =0 (A21) 

5. Example in the 'strong disorder approximation' 

Since the general discussion is rather obscured by the absence of an explicit expression for the eigenvalue Ao(/3, W) 
of the kernel of Eq. A15, it is useful to consider the following strong disorder approximation (called 'upper limit' 
condition in [9]), where the recursion for the X„ in Eq. A5 is simply replaced by [9] 

X„ ^ e„ (A22) 

The argument is that in the limit of very large W, the distribution of X has also a width of order W , so that the 
neglected terms in Eq. A5 are of order K/W. Of course the approximation of Eq. A22 is not very well controlled, 
because it concerns random variables, and it suppresses important correlations between the variables {X, Y) of Eq. 
A5. Nevertheless, it is useful to consider it before returning to the true recursions of Eq. A5, because with Eq. A22, 
the eigenvalue Ao{/3,W) is replaced by simple expression [9] (see also section 7.2 of [17] where the approximation of 
independence between {X, Y) is considered for the Cauchy case) 

Af ^(/3, W) = kJ dep{e)\e\-^^ (A23) 

The velocity then reads 

t;^^(/3,l^) = -^ln(if J dep{e)\e\-^^^ (A24) 
For the flat distribution p(e) of random site energies (Eq. 4), one thus obtains 

The function v^^{f3) is defined on the interval < /3 < 1/2 : it flows towards (-co) in the limits /? — > and /? — > 1/2 
with the following behaviors 

vS''{P,W)^:^^-^ (A26) 

and 

.-(Aiy)^.^^-lln^ (A27) 

It has a single maximum at Pg^^^, given by the extremum condition of Eq. A19, so that /Jf^^.^ is actually independent 
of W. For K = 2, the selected exponent is of order /Jfg^^ — 0.3133 for any W. We note that is value is close to the 
value piyV = 100) c± 0.33 that we measure for the large disorder W = 100 (see Figure 5 a). In conclusion, this 'strong 
disorder approximation' allows to see explicitly how things work on a simple example, and seems to give a reasonable 
value of (iseiec for very large W. 
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6. Argument in favor of < Pseiec{W) < 1/2 



We now return to the analysis of the fuU problem of Eq. A5. Here we should say that we do not agree with the 
discussion of [9] concerning the selection of the tail exponent /3 : 

(i) in [9], the authors conclude that within the localized phase, (3seiec{W) decreases from unity and reaches Pc = 1/2 
exactly at criticality (see the text between Eqs. 6.8 and 6.9 in [9]). 

(ii) in our numerical results of section II we have found instead that the selected exponent Pseiec{W) is always 
smaller than 1/2 : it slightly grows as the disorder strength W decreases and has a value close to 0.5 near criticality 
(see Figure 5 a). 




FIG. 11: Shape of the velocity v{l3, W) as a function of the tail exponent /3 for various disorder strength W : (i) in the localized 
phase W > Wc, the velocity is extremum at a VF-dependent value fHaeiedW) < 1/2 and it is positive Vseiec > (ii) at the 
critical point W — Wc, the velocity is extremum at /3c = 1/2 where it vanishes Vc — (iii) in the delocalized phase, there 
exists another traveling-wave moving in the other direction of velocity Xi{f3,W) = —v{f3,W) > : this velocity is extremum 
at /3 = 1/2 for any W < Wc (see Appendix B for more details) 



We propose the following argument to justify our finding < PseiedW) < 1/2. We think that the tail analysis is 
actually well defined only on the interval < < 1/2 for the following reasons. The eigenvalue Eq. A15 relates the 
behavior of A{X) at |A| ^ oo to the behavior of A{X) near the origin X 0. For A — > cx), the function Q{X + 

will be finite in the region where (A + -^) is finite, i.e. the integration is dominated by the region Ai ^ ^1/A and 
we obtain the power-law decay 

-^(^)|.Poo Ao(^)wU 

However the function A{X) has to be integrable at |A| — > cx) to obtain from Eq. AlO the probability of the variable 
y alone at large y 

f dXA(X) 

P*{y)= dXP*{X,y) ^ ■'-^ (A29) 
According to Eq. A28, the function A{X) is integrable only for 

< /3 < i (A30) 

i.e. we obtain that the tail analysis has a meaning only for < /3 < 1/2. In the limit /3 0, we expect from Eq. A16 
the same behavior as in Eq. A26 

Exactly at (3 — 1/2, the authors of [9] have argued that 

[dMP,W%^,,^^Q (A32) 
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for any W, as a consequence of the symmetry Ao(/3, W) = Ao(l — /?, W) coming from the consideration of the adjoint 
kernel (sec more details around Eqs. 6.7 and 6.8 in [9]). Here in contrast to [9], we think that the region /? > 1/2 is not 
physical because of Eq. A28, but the condition of Eq. A32 is useful to understand why the critical point determined 
by Eqs A21 corresponds to the tail exponent [9] 



(3{W) ^ Jc = l 



(A33) 



To summarize the selection mechanism, the shape of the velocity v(/3, as a function of the tail exponent /3 in 
shown on Fig. 11 for various disorder strength W : 

(i) in the localized phase W > Wc, the velocity derivative is negative at /3 = 1/2 : 



/3=l/2 



^v{P,W>Wc) 



< 



(3=1/2 



(A34) 



The velocity is extremum at a H^-dependent value f^seiedW) < 1/2 and the corresponding selected velocity is positive 

Vselec > 0. 

(ii) at the critical point W = Wc, the velocity is extremum at /3c = 1/2 where it vanishes Vc = 0. 

(iii) in the delocalized phase, the velocity derivative is positive at /3 = 1/2 : 



[d0v{(},W<Wo)]^^,/, 



1 



,2 lnAo{p,W<Wc) 



/3=l/2 



V{0, W < We) 



> 



(A35) 



/3=l/2 



We explain in Appendix B that in this delocalized phase, there exists another traveling-wave moving in the other 
direction with the velocity Ai(/3, W) = — u(/3, W) > : this velocity Xi{l3,W) is then extremum at /? = 1/2 for any 
W < Wc (see Appendix B for more details). 

As a final remark, we believe that the solution exactly at Wc is not valid anymore because of Eq. A28, and that 
the appropriate treatment exactly at criticality should replace the finite velocity motion in vn assumed in Eq. A6 by 
the form 



lny„ = -{cst)n'' + lnt/„ 
where the anomalous exponent < p < 1 has been discussed in the text (see Eqs 35 and 42). 



(A36) 



7. Conclusion 

In conclusion, the analysis of [9] is thus very close to the traveling wave analysis of the directed polymer model 
[26]: the difference is that in the directed polymer, the variables Xn arc random variables independent of the Yn and 
the selected (iseiec can be simply obtained from the non-integer moments of X (see the section 7.2 in [17]), whereas 
in the present localization model, the variables X„ and Yn are correlated and one has thus to solve the eigenvalue 
problem for the integral kernel of Eq. A15. In the following Appendix, we explain how the same ideas can be used in 
the delocalized phase. 

APPENDIX B: TAIL ANALYSIS FOR THE TRAVELING WAVES IN THE DELOCALIZED PHASE 



As explained in section III, in the delocalized phase, we have considered the recurrence for the real Riccati variable 



K 



Xn 



Xn-i{m) 

1=1 ^ ' 

together with the recurrences of Eqs 67 for the auxiliary variables C^^^ introduced in Eq. 66. 



(Bl) 
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1. Tail analysis leading to the same integral kernel as in Appendix A 

To make more visible the similarities with the previous Appendix, it is convenient to perform the change of variables 

Di'') ^ {XlyC^^) (B2) 
to obtain the following recursions for these variables 

The recursion of Eq. Bl for X„ alone will as before converge to some distribution P*{X) satisfying Eq. A7. in the 
delocalized phase, one expects that the variables C*^"^^ will grow exponentially (see Eqs 82). We thus set 

£,(9) =enA,^(g) ^B4) 

where Ag > governs the exponential growth and where d^'^ remains a finite random variable upon the iteration 

m=l ^ \ // 

The similarity with the discussions of Appendix A is now obvious (see Eqs A4 and A5). Assuming the power-law 
decay 

with < < 1, one finds that the function $(X) has to satisfy the eigenvalue problem 

=K j j^Q{X + ^MXi) (B7) 

i.e. it is the same integral kernel as in Eq. A13 with the correspondences /3 — > g/i and e~'"^ — > e^^^", but now we are 
interested in the phase Ag > 0. 

2. Selection of tail exponents 

We have argued above (see section A 6) that the problem for the function $(X) is well defined only for < f3 = 
qfj, < 1/2, i.e. the exponent jjq that governs the power-law of Eq. B6 is restricted to the interval 

< /xg < ^ (B8) 

Denoting as before Aq{P,W) the maximal eigenvalue of the kernel of Eq. A15, each mode /U is associated to the 
Lyapunov exponent (see Eq. A4) 

Xg{f,) = -lnAo{qij,W) (B9) 
M 

Prom Eq. A16, we have the following behavior for — > 

X,{^i) ^^^\nK (BIO) 

Near the other boundary l/{2q), the property of Eq. A32, yields that the partial derivative with respect to fj, is 
still negative at /j. ^ 1/(2(7) 



4lnAo(l/2,l^) 



< (Bll) 

/*^l/(2g) 
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FIG. 12: (a) Probability distribution Hi of the variable lnd^''=^^ of Eq. B4 for various disorder strength W = 1,5, 10, 15, 16.95 
: the exponent /il'^''^'^ of Eq. B6 is of order /if'^''^'^ ~ 0.5 in the whole delocalized phase, (b) For W = 1, comparison of the 
probability distributions Hi and H2 of the variable In d'''' for q = 1 and q = 2 : the exponents of Eq. B6 are respectively 



• 0.5 and /i^" 



0.25 in agreement with Eq. B12. 



For the special case q — 1, the curve Ai(/3) is directly related via \i[f3) — —v{(3) to the curve w(/3) discussed in 
Appendix A and shown on Fig. 11 (see the delocalized case W < Wc). Since the selected tail exponent i^seiec has to 
extremize Xq{fi), we conclude that the selected tail exponents are given by the simple values 




(B12) 



in the whole delocalized phase. This is in agreement with our numerical results. We show on Fig. 12 (a) the probability 
distribution of ff^^i (In d) of the variable (In d) of Eq. B4 for various disorder strength W = 1,5, 10, 15, 16.95 : the tail 
exponent ^f^'-'^'^ of Eq. B6 remains of order ^f^'-'^'^ ^ 0.5. For q = 2, we find similarly that the selected tail exponent 
remains the same within the whole delocalized phase and that it is of order ^^^^'^'^ ^ 0.25. On Fig. 12 (b), we compare 
for W = 1 the probability distributions Hq=ii\nd) and -ffg=2(ln(i) of the variables Ind'^^^^^ : the exponents of Eq. 
B6 are respectively ~ 0.5 and ^2^^"^^ ^ 0.25. 

From Eq. B12, we obtains that the Lyapunov exponents satisfy the identities 

Xf'"' = qXf'"' (B13) 

As explained in the text around Eq. 85, this identity is important to understand the decay of I.P.R. with the number 
N of generations in the delocalized phase. 



3. Finite-size corrections introduced by the pool method 



As explained in section II C 1, the pool method consists in representing a probability distribution P{x) by a large 
number Mpooi of random variables {xi}. In particular this introduces a cut-off in the tail of P{x) around Xmax with 
P{xjyiax) ~ ^/Mpooi- In the field of traveling waves, the presence of such a cut-off has been much studied (see [43] 
and references therein) : the leading correction to the selected velocity is usually logarithmic in Mpooi 

Vselec{Mpool) - Vselec{oo) ^ — (B14) 

(in Mpooi) 

For our present problem, where the critical point corresponds to a vanishing velocity w = 0, we thus expect that the 
true critical point Wc{oo) corresponding to Vseiecioo) = will be shifted towards a pool-dependent pseudo-critical 
point Wc{Mpooi) corresponding to Vseiec{Mpooi) = 0. As explained in the text, all the traveling waves encountered in 
the present paper vanish linearly in {W — Wc). We thus expect that the difference [Wc{Mpooi) — Wc{oo)] will also 
decay only logarithmically as 1 / {In Mpooi)" ■ From the analysis presented in Appendix A and B, it is obvious that 
the pool dependent critical point for the traveling wave of the Landauer transmission (see Appendix A) and for the 
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auxiliary variable C^^"^ or D^^^^ (sec Appendix B) arc the same because they are defined by the same condition in 
terms of the eigenvalue Aq{0) of the integral kernel. However, the other variables C*^*^ or D^'^^ that have the same true 
critical point Wc(oo), will not have the same pool-dependent critical point Wc{Mpooi) , because the constant J in Eq. 
B14 will depend on q. This is indeed what we observe with our numerical computations with the pool Mpooi = 10^ : 
the pseudo-critical point Wc{Mpooi) for the variable C^^^ is below the pseudo-critical point VFc(-Mpooj) for the variable 
C^^\ We observe that the difference between the two becomes smaller for the bigger pool Mpooi = 10^. In conclusion, 
numerical studies based on the pool method are valid for observables that are related to a single traveling wave, but 
one cannot study the critical properties of observables that depend on two distinct traveling waves that have different 
pseudo-critical points. This is why in the text, we have not been able to present reliable numerical results for the 
statistics of I.P.R. I2 that depends on both C^^^ and C^^^. 



APPENDIX C: REMINDER ON THE TRAVELING/NON-TRAVELING PHASE TRANSITION 

STUDIED IN REF. [30] 

In this Appendix, we briefly recall the finite-size scaling properties of the traveling/non-traveling phase transition 
studied in Ref. [30], because these properties are useful to interpret our numerical results described in the text. 

For a one-dimensional branching random walk in the presence of an absorbing wall moving at a constant velocity 
V, the survival probability Q{x, t) presents a phase transition a,t v — Vy with the following critical behaviors (see more 
details in [30]) : 

(i) for V < Vc, it converges exponentially in time towards a finite limit Q*{x) > 

Q{x,t)^Q*{x)+e-icl){x) (CI) 
where the limit Q* {x) presents an essential singularity 

Q*{x) ~ e~T^^f^ with k=^ (C2) 
and where the relaxation time r diverges as 

1 3 

r oc ^ with v = - (C3) 

(ii) for V > Vc, it converges exponentially in time towards zero 

0(x,t)~e-* (C4) 

where the relaxation time f diverges as 

f oc — with i/ = 1 (C5) 

(iii) exactly at criticality, it converges towards zero with a stretched exponential 

Q(a;,i) -e-^'^*^)*' with p=\ (C6) 

(iv) the critical region is described by some finite-size scaling form governed by the relaxation time r that appears 
in Eq. CI (and not by the relaxation time f appearing in Eq. C4) 

\nQ{x,t)ci-tPG[t^^''{vc-v)^ (C7) 

where the scaling function G{u) has the following asymptotic behaviors. For u — » +00, has the one has the power-law 

G{u) ~ u-'P (C8) 

to recover the finite limit of Eq. C2 and one has the scaling relation 

K = iyp (C9) 
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For u — > — 00, one has the power-law 

G{u) ~ {-uY^^-P^ (CIO) 

u—* — oo 

to recover the exponential decay of Eq. C4, and one has the scaling relation 

i> = u{l-p) (Cll) 

For the problem considered in [30] all the exponents are exactly known p = 1, p = 3/2, p = 1/3 and k = 1/2. 

In our numerical results presented in the text concerning the traveling/non-traveling phase transitions that occur 
in Anderson localization on the Bcthc lattice, wo find very similar critical behaviors : the velocity of the traveling 
wave vanishes linearly with p = 1, and the finite-size scaling in the critical region is governed by the other exponent 
v that appears in the relaxation towards the finite value of the non-traveling phase. 
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